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ABSTRACT 



Aims. We investigate the effects of radiative shocks on the observed X-ray emission from the Type II supernova SN1993J. 
To this end, the X-ray emission is modeled as a result of the interaction between the supernova ejecta and a dense 
circumstellar medium at an age of 8 years. 

Methods. The circumstances under which the reverse shock is radiative are discussed and the observed X-ray emission 
is analyzed using the numerical code described in Nymark et al (2006). We argue that the original analysis of the 
X-ray observations suffered from the lack of self-consistent models for cooling shocks with high density and velocity, 
leading to questionable conclusions about the temperatures and elemental abundances. We reanalyze the spectra with 
our numerical model, and discuss the expected spectra for different explosion models for the progenitors. 
Results. We find that the spectra of SN 1993J are compatible with a CNO-enriched composition and that the X-ray 
flux is dominated by the reverse shock. 
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1. Introduction 

Many core-collapse supernovae (SNe) show strong X-ray 
and radio emission. This is, in particu lar, true for the 
Type IIL, Iln and lib SNe (see e.g. Ilmmler &: LewinI 
120031 for a review). We now know that this is a re- 
sult of the interaction of the supernova ejecta and the 
circumstellar medium (CSM) of the supernova. A two- 
shock structure is formed consisting of the outgoing su- 
pernova shock (hereafter the forward shock) an d a reverse 
shock, which moves backward into the elect a (IChevalieJ 
fT982allbl : [Chevalier fc FramsOTlfT99l hereafter [CFg^Trhe 
region behind the forward shock consists of shocked cir- 
cumstellar material, with T ~ 10* — 10^ K, while the re- 
gion behind the reverse shock, which consists of shocked 
ejecta, is cooler, with T ~ lO'' K. The radio emission is 
cause d by synchrotron emiss ion from the forward shock 
(e.g. iSramek fc Wieil [200l . while the X -ray emission 
originates partly from free-free and inverse Compton emis- 
sion behind the forward shock, and partly from line emis- 
sion in the cooler gas behind the reverse shock. Both the 
radio and the X-ray emission are powerful tools for under- 
standing the physics of both the SNe and the CSM. 

With the new generation of X-ray telescopes, XMM- 
Newton and Chandra, more detailed spectra have been 
obtained than ever before. This puts greater demands on 
the tools used for analyzing the spectra. Traditionally, X- 
ray spectra have been analyzed with the XSPEC package, 
which includes a n umber of models that can be co mbined to 
get the best fit. In lNvmark et al.l (|2006l . hereafter |NFK06D 
we showed, however, that the SNe having the strongest X- 



ray emission are likely to have radiative shocks. This has 
not been accounted for in the analysis of the SNe observed 
in X-rays. Furthermore, none of the models included in the 
XSPEC distribution offers the possibility of fitting cooling 
shocks in a consistent manner. 

In radiative shocks a wide range of temperatures con- 
tribute to the spectrum. The cooling also affects the hy- 
drodynamic structure and influences the emitting volume, 
and thus the total luminosity from the interaction region. 
This is not properly accounted for when fitting with one, 
or even several, single-temperature models, and the el- 
emental abundances can be over - or under-estimated in 
the resulting analysis. In iNFKOGi we described a numeri- 
cal model that computes the emission from a cooling shock 
self-consistently. This was applied to t he coronal emission 
from t he ring collision in SN 1987A by iGrdningsson et al.l 
(|2006D . Here we apply this model to the X-ray spect ra of 



SN1993J, obtai ned with Chandra (ISwartz et ai.ll2003f) and 
XMM-Newton (jZimmermann fc Aschenbachll2003[ r 
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2. Observations and data analysis 

SN 1993J was discove red on March 28, 1993 in M82 
(jRipero fc Garcial I1993D . It was one of the brightest su- 
pernovae of the twentieth century, and has been followed 
in great detail in many wavelengths in the decade fol- 
lowing the explosion. On the basis of the detection of 
hydrogen lines in the spectrum SN1993J was first iden- 
tified as a type II supernova, but as the Balmer lines 
soon weakened its spectrum came to resemble that of a 
type lb supernova, and it was re-classified as a Type lib 
(|Filippenko et al.]|l993t) . This transition is explained by an 
almost complete loss of the hydrogen envelope due to heavy 
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pre-supern ova mass loss, probably due to a binary interac- 
tion (e.g . 'Nomoto et al."1993'; 'Podsiadlo wski et all Il993t 
iRav et al. 1993; Woo slcy e t al. 19 94). This hypothesis was 
confirmed bv lMaund et al.i (j2004f ). when they detected the 
unambiguous signature of a massive star, i.e., the binary 
companion of the supernova progenitor, in HST observa- 
tions of the supernova 10 years after the explosion. Evidence 
of circumstellar intera ction was seen at early times in radi o 
(IVan Dvk et al.lll994D. UV (iFransson fc Sonnebornl ll994D. 



X-rays ( Zimmerman et all ll994t"Uno et al.l 20021 ) and 7- 
rays by OSSE on CGRO (iLeising et all ll994DrThe detec- 
tion of the supernova in the radio (IChandra et al" l2004t 
IWeiler et al ] i2007, ) and X-rays bands (|Uno et al.ll2'00g ) even 
after a decade shows the presence of an extended dense 
GSM. 

During the first weeks ROSAT and ASCA observations 
showed a hard X-ray spectrum. This was confirmed by 
OSSE, which showed a hard spectru m, well fit by a free - 
free spectrum with kT - 100 keV (|Leising et al.lll994l ). 
However, after six months ASCA observations showed that 
this had changed to a soft spectrum (Uno ct al. 2002). This 
transition is consistent with circumstellar interaction, as 
the decreasing optical depth in the expandi ng cool shell al- 
lows m ore of th e soft em ission to escape (jFransson et al.l 
(fT99l . hereafter [FLCH). 

In May 2000 S N 1993J was observed with Chandra 
(jSwartz et al.ll2003[ ). Strong Fe L emission is evident in the 
spec trum, as are emissio n lines from Mg XI-XII and Si Xiii- 
XIV. ISwartz et al.l (|2003l ) fit the spectrum with three com- 
ponents - two low temperature components (kT ^ 0.35 keV 
and kT ~ 1.01 keV) with no He, C, O and Ne but enhanced 
in N, Mg, Si and Fe; as well as a high temperature com- 
ponent (kT ^ 6.0 k eV), consistent with solar abundances 
ijSwartz et al]|2003l ). This was argued to be consistent with 
an interaction model, in which the high temperature com- 
ponent comes from the circumstellar shock, while the two 
low-temperature components come from different parts of 
a radiative revers e shock. 

In April 2001 IZimmermann fc AschenbachI (|2003D ob- 
served SN 1993J with XMM-Newton. Their spectrum was 
found to be best fit with two components (with kT ~ 
0.34 keV and kT ^ 6.4 keV). These are in good agree- 
ment with the high and low temperature components 
of the Chandr a spectrum, b ut th e intermediate compo- 
nent found by ISwartz et al.l (j2003D in the Chandra spec- 
trum is missing in the analy s is of the XMM spectrum by 
IZimmermann fc AschenbachI ()2003l ). This could be caused 
by a real change in the spectrum between the two observa- 
tions. Also here the two components are taken as evidence 
of a two-shock structure, but they found the abundances 
to be poorly constrained, and no conclusions were drawn 
about the composition of the ejecta. 

As we argue below, both these analyses are doubtful if 
the shock is radiative, as it most likely is. The abundances 
were kept as completely free parameters, and it is not clear 
how unique the solution is, especially with a realistic tem- 
perature model. The thin hydrogen envelope means that 
the reverse shock will traverse the hydrogen rich region on 
a short timescale, and the emission will already after a year 
be dominated by processed material, which gives us an ex- 
cellent opportunity for testing stellar evolution theories. In 
this paper we take a different approach and use abundance 
models based on self-consistent stellar evolution models as 
input. 



For the analysis in this paper we make new extractions 
of both the Chandra and XMM-Newton data to fit with 
our models. Here we describe the observational and data 
analysis procedure. 

2.1. The Chandra dataset 

The parent galaxy of SN 1993J, M81, was observed with a 
50ks exposure on May 7, 2000 with the Chandra Advanced 
CCD Imaging Spectrometer (ACIS) spectroscopy array op- 
erating in imaging mode [Swa rtz et al. 2003). We refer to 
ISwartz et al.l ()2003f ) for the technical details of the observa- 
tions. We extracted the data from the ChandraXO archive 
and analyze d this using s tandard CIAO analysis threads 
and XSPEC l|Arnaudl ll996). In order to check for any X-ray 
flaring events, we selected a background region of size 90", 
carefully excluding point sources. We used a radius of 7" 
centered at the SN 1993J position to extract the SN events; 
and an annulus of inner radius 7" and outer radius 18" cen- 
tered at the super nova position to obt ain t he background 
count s. We refer to ISwartz et al.l ()2003D and lChandra et al] 
(|2005D for the data analysis procedure. We binned the data 
in 15 cts per channel for the spectral analysis. The observed 
spectrum is shown in Fig[T] While analysing the supernova 
spectrum with XSPEC, we ignored the counts below 0.2 
keV and above 6.0 keV, since there was almost no flux in 
these energy ranges. 

2.2. The XMM dataset 

XMM-Newton observed SN1993J on April 22/23, 2001, for 
a total of 132 ksec. The 90 ks data was with the PN camera 
in the small window mode and 83 ks with the M0S2 camera 
in image mode. The medium filter was applied both in the 
PN and the M0S2 instrument setup during the observation 
(jZimmermann fc Aschenbachll2003f ) . Due to the presence of 
the variable particle background in parts of the observing 
period, we flagged the flaring events from the light curve. 
In the PN observations about half of the observations were 
affected by bursts in the particle background. We extracted 
the spectrum and created response matrices using XMM- 
SAS. The spectrum was binned as group of 50 bins. The 
observed spectrum is shown in Fig. [TJ Spectral fits were 
performed using only the PN camera data offering the best 
statistics. 



3. Modeling 

In general, the X-ray flux is a combination of the flux from 
the forward and reverse shocks. As we will show below, in 
most cases the reverse shock is likely to be radiative. For 
these cases we have made a self-consistent model in which 
hydrodynamic calculations for a stationary flow are com- 
bined with time-dependent ionization bala nce and multi- 
level atoms, including optical depth effects l|NFK06l ). This 
means that the emission in the cooling flow behind the re- 
verse shock is followed from the shock to the cool shell be- 
hind the shock, giving correct weight to the emission from 
each zone. The parameters of the shock calculations are 
the reverse shock temperature, the density and the chem- 
ical composition, and the output is the emitted spectrum. 
For the cases where the cooling time is long we assume a 
one-temperature spectrum. 
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Fig. 1. The Chandra dataset (left) at t=2594 days since explosion and the XMM dataset (right) at t=2943 days since 
explosion. 



Ionization rates are from Arnaud fc Rothenflud (jl985[ ) 
and lArnaud fc Raymond ( 1992), while recombination rates 
arefromlNahan(il996.il999D.lNahar fc PradhanI (|1997D. and 
IZatsarinnv et al.l (|2003l l2004allbl ). Most of the collisional 
and radiative data are from the Chianti database, version 
4.2 (|Dere et al.l 119971: lYoung et al.ll2003f FI. For further de- 
tails we refer to lNFKOd 



Houck & Fr anssonI 1199^ 
the 4H47 model 



M 



4. Structure of the ejecta 

The properties of the circumstellar interaction depend on 
the density structm'e in the ejecta. Explosion models for 
red supergiants indicate that the density distribution in 
the ejecta can be described by a power law, p oc . For 
a polytrope of ind e x n = 3/2, typical of red supergiants, 
iMatzner fc McKed ([1999') have shown that 77 « 10 in the 
outer layers. Deviations from a simple polytropic progeni- 
tor may, however, result in a more complicated structure. 
As an exam ple, the model 4H47, specifically calculated for 
SN 1993J (jShigevama et al.l 119941 : iNomoto fc Hashimotol 
Il988[) , shows a structure with a steep outer part and a flat- 
ter inner profile, and then an additional steepening further 
in, as well as considerable structure in the inner parts (see 
their Fig. 3). This profile is , however, also a simplification, 
since llwamoto et al.l ()1997f) find that especially the inner, 
steep region is Rayleigh- Taylor unstable. 

As the reverse shock moves inward in the ejecta, it 
passes slowly through the steep outer gradient, while mov- 
ing more rapidly in the flatter part before eventually reach- 
ing the second steep gradient. When the gradient changes, 
the spectrum is also expected to change as the reverse 
shock slows down (steep part) or speeds up (flat part). 
There is evidence that this has happened in SN 1993J 
(|Dwarkadas et al.ll2005f) . 

The shock radius depends on time as Rs oc t™, where _2 
the expansion parameter m = (p — 3)/ ir] — 2). Numerical -^'^rcv = 1.3 x 10 (77 — 4) 
calculations bv lDwarkadas et all ()2005f ) based on the 4H47 
model indicate that a change in the expansion parameter 
should occur at approximately day 2300. Such a change 



was indeed observed by iBartel et al.l (|2002f ) in the radio 
at ^ 1400 days. For the cases discussed in the present 
work it is worth noting the differenc e in the X-ray spec- 
tra from Chandr a (ISwartz et al. 2003f) on dav 2594 and the 
one from XMM (Zimmcr mann fc Aschenbachll2003[ ) a year 
later. Although both these observations were done some- 
what later than the expected change in the spectrum, the 
observed difference could be related to a change in the ex- 
pansion parameter. 

Models of the light curve and spectrum of SN 1993J 
strongly indicate that less than 0.5 Mq of the hydro- 
gen envelope was left at the time of the explosion (e.g . 



Wooslev et al.lll994 iNomoto et al.lll993l : iRav et all [19931 : 



Further, explosion model s like 
Nomoto fc Hashimotol ()1988D and 



IShigevama et al.l (|l994f ) predict that the outer parts of the 
ejecta should contain nearly equal mass fractions of H and 
He as a result of CNO processing. 

Assuming a constant density gradient in the ejecta and 
a steady wind, we can estimate the swept up mass and 
therefore what the composition at the reverse shock should 
be. If the density of the CSM is given by 



the mass swept up by the reverse shock is 

_ r/ - 4 M 
2 



(1) 



(2) 



where is the wind velocity ()Chevalierj|1982b| ). With Rs — 

Kji, we get 



104 km s- 



(3) 



4 / V 1000 days 



Mq, 



^ CHIANTI is a collaborative project involving the NRL 
(USA), RAL (UK), and the following universities: College 



where Vc\ is the 



London (UK), Cambridge (UK), George Mason (USA), 
Florence (Italy). 



ejecta velocity and the time 
and smce explosion in days. For SN 1993J = 
(M/10-5 Mq yr-i)/(uw/10 km s'^) w 4 (|FLC9fiD . 
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For the Chandra observations we take = 2600 and 
V4 = 14j/10^km s~^ « 1, which is consistent with the 
expa nsion velocity derived fro m VLBI observations at this 
time (jBartel et al.ll2OO2ll2OO70 . 

Assuming an ejecta density gradient 77 = 25 in the outer 
zones, and 77 = 5 in the shallow middle zone, as indicated 
by the 4H47 model by Shig evama et alj (|1994f ). we find 
the swept up mass at 2600 days to be ~ 0.38 Mq. The 
hydrogen-rich envelope in the 4H47 model has a mass of 
0.47 Mq. Inside this envelope the outer part of the He zone 
is N-rich, while the inner part is C-rich. The He/N zone 
contains 0.3 Mq, and the He/C zone contains 1.5 Mq. 

The estimate above shows that the reverse shock may 
be close to the border between the H-rich envelope and the 
He-rich zone. However, this estimate depends on the mass 
assumed for the hydrogen envelope . With a lower hydro - 
gen mass, as m the sl5s7b model bv lWooslev et al.l (| 19941) . 
where the hydrogen envelope is only 0.2 Mq, it is therefore 
quite possible that the reverse shock is somewhere in the 
He-rich zone. Nevertheless, we cannot exclude the possibil- 
ity that clumps of metal-rich material may have reached the 
outer parts by mixin g, as indicated by the simulations by 
llwamoto et al.l ()1997[ ). In Sect. [5] we therefore also discuss 
th e spectrum resulting from metal-rich compositions. 

ISwartz e t al.l 1)200 3') found that the Chandra spectrum 
of SN 1993J was best fit with a composition dominated 
by H, and where N, Mg, Si, and Fe were enhanced with 
respect to solar abundances, while He, C, O and Ne were 
absent. These abundances are, however, not consistent with 
the calculated models for t he ejecta. In the analysis by 
IZimmermann fc AschenbachI (|2003) the XMM observations 
were fit with a number of different models, all with widely 
differing elemental abundances, but these abundances were 
very poorly constrained. 



5. Radiative shocks 



The structure of the reverse shock was discussed in lNFK06l . 
where we showed that this can remain radiative for several 
years, and longer. The emission and cooling depend on the 
reverse shock temperature, the pre-shock density and the 
composition. These are the fundamental parameters, but 
each of these are related to the shock velocity, Vg, mass 
loss parameter, A,, ejecta density gradient, 77, and the time 
since explosion, t, by the relations discussed in NFK06. For 
clarity we here repeat the relevant relations. 

An estimate of whether the shock is radiative can be 
found by co mparing t he expansion time with the cooling 
timescale. In iNFKOGi we showed that the ratio between 
these timescales can be written as 



^cool 

t 



6.7 X 10-2^1^^4 id 



ifiA 2)2(7] - 3) (77 - 4)A(Trcv)i* ' 



(4) 



where /i is the mean mass per particle and /ia is the mean 
atomic weight. The cooling function, A(Ti.ov) is discussed in 
INFK06 . The shock is radiative if icooi ^ t. If the interaction 
region is thin compared with the size of the supernova, 
as will be the ca se if the shock is radiative, the thin shell 
approximation of lChevalieil ()198 2b) applies. In this case the 
temperature of the reverse shock is given by 



Using this in Eq. ^ and requiring that tcooi/t ^ 1.0 we find 
that the shock is radiative if the reverse shock temperature 
fulfills the requirement 



X 10 
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Ma(?7 - 2) 

(/z(/.A - f^)iv - 3)(7? - 4)A(T)i,tji 



(6) 
(7) 



_ 2.27X 10V ^^2 
(r;-2)2 



(5) 



Table [2] shows the limiting temperatures, Tjad, and ex- 
pected shock temperatures for a number of compositions 
and density gradients for the case of SN 1993 J, with A^, — A 
and id — 2600 days. T^ad is defined as the limiting temper- 
ature corresponding to icooi = t, such that if T^cv < Trad 
the shock will be radiative while Trev > Tiad gives an adi- 
abatic shock. These regimes are also illustrated in Fig. [5] 
in terms of the ejecta density gradient and the shock ve- 
locity. The lines in Fig. [2] represent the border between the 
adiabatic and radiative regimes for a few different composi- 
tions. If the ejecta density gradient and the shock velocity 
are in the region above the line corresponding to the rele- 
vant composition, the shock is radiative. This means that 
for a given shock velocity, the shock is more likely to be 
radiative the higher the density gradient is. We emphasize 
that the values of 77 and Vg corresponding to radiative and 
adiabatic shocks shown in Fig. [2] are a function of the time 
since the explosion, the mass loss rate and the cooling rate, 
and are therefore only valid for the specific case discussed 
here^ 

ISwartz et all (|2003 D fit the SN 1993J data with two low 
temperature components (fcT = 0.35 keV and kT =1.0 
keV) and one high temperature component {kT — 6.0 
keV). From Table [2] we find that the fcT^v = 0.35 keV 
shock is likely to be radiative, and, depending on the com- 
position, also the fcTrev = 1-0 keV component, while a 
fcTrev = 6.0 kcV schock is probably adiabatic. This means 
that the assumption of a single-temperature shock is not re- 
alistic for the two low-temperature components, and their 
analysis is therefore not consistent. 

It is possible that the high energy component comes 
from the hot outer shock, while the fcTrcv = 1.0 keV compo- 
nent comes from a radiative reverse shock. In this case, the 
kTj-cv = 0.35 keV component could be representative of the 
cooling part downstream from the reverse shock. However, 
as we discuss below, in Sect.[6l the temperature of the outer 
shock is unlikely to be as low as 6.0 keV. Another possibilty 
is the presence of several shocks in the ejecta, for example 
if the ejecta is clumpy. In that case all three components 
could come from the reverse shock, but corresponding to 
different angles with respect to the clumpy ejecta. This al- 
ternative is discussed in more detail in Sect. [71 

6. Emission from the circumstellar shock 

At late times the reverse shock is expected to dominate 
the X-ray emission. There may, however, still be a con- 
tribution from the outer shock, in particular at energies 
^ 5 keV. The high temperature component found in the fits 
to the Chandra and XMM observations of SN 1993J was 
attributed to the circumstellar shock, which is expected to 
have a considerably higher temperature than the reverse 
shock. It is surprising, though, that this temperature is as 
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lo w as ~ 6. keV. The shocked CSM has an ion temperature 
of (|FLC96D 



2.27 X IOVt^ 

vn 



if 



Vf K, 



(8) 



where = 0.61 for a solar composition, and given in 
Table [T] for other compositions. For V4 = 1 and 77 = 12, 
kTca <: 100 keV. Allowing for the possibility of a lower 
ejecta density gradient or shock velocity, this value may 
be reduced somewhat, but it is difficult to get a tempera- 
ture lower than ~ 40 keV if the ion temperature is equal to 
the electron temperature. 

However, if coUisionless heating is inefficient, the elec- 
tron temperature can be considerably lower than the ion 
temperature (e.g., lFLC96f ). A lower limit to the electron 
temperature may be found from the timescale for electron- 
ion equilibration 



An electron t empe rature as low as 6.0 keV, as found by 
ISwartz et al.l ()2003D . is therefore in principle possible. Also 
cosmic ray dominated shocks can give a lower temperature 
and higher compre ssion rate (e.g.. lDecourchelle et al.|[2000l : 
lEUison et al.ll2004l ) , which could contribute to lowering the 
temperature of the forward shock. For efficient cosmic ray 
acceleration, with an injection efficiency (fraction of protons 
with superthe rmal energies) of 77 = 10~- ^ at the outer shock 
the models of iDecourchelle et al.l (|2000D give temperatures 
of a few keV, while the resulting high compression can lead 
to relative thickness of the outer shock (as defined below) 
of / « 0.04. Also this effect can therefore give a shock 
temperature of 6.0 keV. 

However, a second requirement is that this compo- 
nent should also produce the observed luminosity of 1.6 x 
Iq3& Qpgg g-i f^iiQ range 2.4-8.0 keV. The luminosity from 
the circumstellar shock is given by 



^e — i 



29 
0.21 



10^ K 



3/2 



10^ cm^3 



V 109 K 



3/2 



A* ^Vit\ days 



(9) 



CoUisionless heating by plasma instabilities may, how- 
ever, increase the electron temperature. There is, however, 
evidence from observations of galactic supernova remnants 
that coUisionles s heating is only effec t ive for low Mac h 
number shocks (jGhavamian et al.l[2007t IVink et al.ll2003^ . 
Because of the much higher densities at these comparatively 
early stages, it is, however, not clear how these results apply 
to our cases. 

Equating tc-i to the expansion time gives 



kT, Z 7AV; 




-2/3 



keV. (10) 



L 



A.cs 



2.54 X 10^^ l^^J^Alvr' ( 1- 

MMA * ^ VO-2 

^y-0.236;^-2g-1.44/AT« 



V 1000 days 



erg s 



(11) 



using a Gaunt factor gs ~ 1.87Tj 



-0.264 



(|FLC96i) . In this 



expression / is the relative thickness of the adiabatic outer 
shock, defined by / = [Res — Red) /Red, where Res is the 
outer shock radius and Red is the radius of the c ontact dis- 
continuity. For 77 = 7 - 20, / « 0.2 - 0.3 (|CF94f ). Inserting 
the appropriate parameters, we find that the expected lu- 
minosity in the range 2.4-8.0 keV from a circumstellar com- 
ponent with kT ~ 6.0 keV is Les ~ 2.0 x 10^^ ergs s"^ at 
2600 days, i.e., only about one tenth of the observed lumi- 
nosity. A higher compression, as resulting from cosmic ray 
acceleration, can lead to higher luminosity. However, even 
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the most efficient model in iDecourchelle et al.l ()200Cll ). with 
77 = 10~^, which gives a compression of / « 0.04, will for a 
temperature of 6.0 keV lead to an increase of the luminosity 
by only a factor of five, which is not enough to explain the 
observed luminosity. We conclude that it is unlikely that 
the circumstellar shock is responsib le for th e 6.0 keV com- 
ponent. This is in agreement with 'F LC96I who find that 
the luminosity of the forward shock after t ~ 1000 days is 
orders of magnitude lower than that of the reverse shock. 



7. Oblique shocks 

If either the CSM or the ejecta is clumpy, parts of the re- 
verse shock will hit these clumps obliquely. The transmitted 
shocks will then have a lower velocity than the part of the 
reverse shock which hits the clumps head-on, and also the 
shock temperature will be lower. This scenario could result 
in a considerable range of shock velocities contributing to 
the emission. In particular, even if the head-on shock is adi- 
abatic for the parameters chosen (i.e., a given reverse shock 
velocity and ejecta density), the tangential shocks may be 
radiative. 

This is similar to the situation found for t he ring col- 
Usion in SN 1987A, where IZhekov et all (|2006D find that a 
range of shock velocities are necessary to explain the ob- 
served spectrum. The presence of both hard X-rays and 
optical emission from the shocked CSM clearly shows the 
importance of geometric and complex dyna mics of a clumpy 
CSM. The hydrod ynamical simula t ions b y llwamoto et al.l 
(|1997f) . as well as iKifonidis et al.l ()2006[ ). however, show 
that also the ejecta is subject to Rayleigh- Taylor insta- 
bilities in connection to the passage of the shock. The 
same situation as in the CSM case with a range of veloc- 
ities for the reverse shock will then arise. With this back- 
ground, we therefore attempt to fit the spectra with two or 
more shocks, possibly with a mix of adiabatic and radiative 
shocks. The results of this fitting are discussed in the next 
section. 



8. Model fitting of the X-ray spectra 

In this section we discuss the modeling of the spectra along 
the lines of the previous discussion. First we emphasize that 
we do not aim at perfect statistical fits, keeping all param- 
eters, like temperature, abundances and absorbing column 
density, free. Instead, our aim is to investigate how well - or 
how badly - the observations are reproduced by different 
realistic explosion models. Therefore we focus on the im- 
portant spectral features, and investigate how well each of 
them is reproduced by each model. In this way we hope to 
distinguish between different models for the composition, 
and therefore the nature of the progenitor. We believe this 
to be more physically illuminating than the more commonly 
used method of varying all parameters until a best fit has 
been achieved. In particular, with data of low spectral res- 
olution and moderate S /N this can easily result in doubtful 
abundances. Different sets of abundances may result in sim- 
ilar values of , and the one with the lowest value may have 
little significance compared to others with substantially dif- 
ferent abundances. Since we do not allow all parameters to 
vary freely, the statistic is not really appropriate. One 
approach which is used in similar circumstances is to inves- 
tigate quantities which can be determined unambiguously 



in the data and the model fe.g.. [Badenes et al.ll2008l ). In 
the present case, however, the quality of the data does not 
allow any this to be determined unambiguously. Therefore 
we choose instead to focus on specific spectral features, and 
use the statistic as a guideline to distinguish between a 
fe w specific models . This is similar to the approach taken 
bv lRakowski et all (|2006[ ). 

We have chosen to relate our models to the ejecta com- 
position foun d from two explosion models for SN 1993 J, the 
4H47 model (jNomoto fc IIa shimoto"19 88l: IShigevama et"al] 
I1994D and the sl5s7b model (Wooslcv et al.lll994f) . From 
each of these we have selected a number of compositions, 
corresponding to the different burning zones in the ejecta 
(Table [1]). 

We also discuss a CNO-enriched composition resembling 
the one in the circumstellar ring of SN 1987A, but with 
initially three times higher metallicity, to account for the 
low metallicity of the LMC compared to the more normal 
metallicity in M81 and NGC 3877. This composition has 
been determined in the following way: We took H/He=0.25 
by number, as appropriate for the circumstellar ring of 
SN 1987A. The sum of C, N and O abundances was taken 
to be equal to the solar value, while the ratios were the 
same as in the ring of SN 198 7A, i.e., N/C=5 and N/0=1.1 
(jLundgvist fc Franssonlll99^ . For all other elements we use 
solar abundances. We note that the optical observations o f 
SN 1993 J indicate a high N/C ratio (jFransson et al.ll2005[ ). 

We also note that due to the presence of the cooling 
shell, part of the fiux d ensity from the reverse shock will 
be absorbed (e.g. ICF94l ). This absorption will come in ad- 
dition to the Galactic absorption column density. However, 
unlike the Galactic absorption, which approximately should 
have solar abundances, this extra absorption should nor- 
mally have the abundances of the burning zone close to the 
shock. We therefore first determine the column density for a 
nearby source in the field of view of the SN, and freeze the 
Galactic column density to that value in our fitting. The 
extra absorption contribution will then solely come from 
the cooling shock. 

For each composition we attempt to fit the spectrum 
with one reverse shock temperature, which we allow to vary 
between the maximum reverse shock temperatures corre- 
sponding to ejecta density gradients in the range = 7 to 
20. The corresponding reverse shock temperature, Tiov, for 
each composition is given in Table [21 together with the lim- 
iting temperature for a radiative shock, Trad- If the best fit 
reverse shock temperature exceeds Tj-ad, we take an adia- 
batic component with a temperature that is allowed to vary 
between Tj-ad and Ticv 

We also tested adding a hot, single-temperature compo- 
nent with solar abundances, representing the outer shock. 
For this component we used the XSPEC model "Vmekal" , 
which gives the spectrum from a hot, diffuse gas with a sin- 
gle temperature. This is basically a free-free spectrum for 
kTcs ^ 10 keV. 

9. Results 

The XMM-Newton spectrum of SN 1993J (Fig. [H right 
hand panel) is dominated by the Fe L emission at 0.7- 
1.0 keV, which is blended with strong Neix-X emission. A 
lower peak is evident at 0.4-0.5 keV, most likely caused by 
strong Nvii lines. Several emission lines are present above 
~ 1.0 keV, in particular from Mgxi-Xii and Sixiii-xiv, 
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with the Mg xi and Si xiii hnes being the stronger, as well as 
Sxvi. The absence of O emission at 0.6-0.8 keV is striking, 
as the O viii lines at 0.65 and 0.775 keV are expected to be 
strong over a wide range of temperatures. 

The Chandra spectrum (Fig. [U left hand panel) shows 
much the same features as the XMM spectrum, but with 
a more pronounced double appearance of the Fe L peak, 
with one peak at ~ 0.8 keV and a second peak at ^ 0.9 
keV. Because the S/N is considerably lower for the Chandra 
spectrum, we will in the following concentrate on the XMM 
spectrum. 

We will now systematically investigate to what extent 
the abundance models in Table [1] reproduce these observa- 
tions. We start with a solar composition, and then go on to 
more specialized models for different nucle ar burning zones. 

In the case of SN 1993 J ISwartz et all ((2003) found the 
best fi t column density for t he M8 1 nucleus to be 4.4 x 10^" 
cm~^. iDickev &: Lockmanl ()1990f ) found that the Galactic 
absorption column density in that direction is 4.17 x 10^° 
cm~^, which is well in agreement with that of the best fit 
column density. Hence we froze the Galactic column density 
with solar composition to this value. 

When fitting the spectrum with one radiative shock, 
we found that the best fit temperature exceeded Trad in 
all models. An adiabatic shock with a higher temperature, 
gave a better fit, but in most cases the best fit was found for 
a combination of an adiabatic shock and a radiative shock 
at half the temperature of the adiabatic shock. 

Adding a hot component from the circumstellar shock 
did not improve the fits, and the best fit temperature for 
this component was ^ 5 keV in all cases. As discussed in 
Sect. [HI the forward shock is unlikely to have such a low tem- 
perature. We therefore conclude that there is no indication 
of a contribution from the forward shock to the spectrum 
of SN 1993J. The absorption column density needed is in 
all cases just a little bit higher than the Galactic column 
density, which suggests that SN 1993J has significant con- 
tribution from the adiabatic reverse shock too, along with 
the radiative part. 

We will now discuss the specific abundance zones one 
by one. The best fit parameters for the solar composition, 
the CNO-enriched co mpositions from the model sl5s7b2f 
(|Wooslev et al.lll994D and the SN 1987A-Hke composition 
are given in Table [31 The spectra were found to be best fit 
with one adiabatic component (denoted here by Ad. Comp) 
and one radiative component (Rad. Comp). For each com- 
ponent the table shows the best fit temperature, and the 
flux normalization constant (Norm) required to fit the spec- 
trtun. 

9.1. A solar composition 

First we note that the relative strengths of the Mgxi and 
Mgxii lines, and of the Sixiii and Sixiv lines, indicate 
that there is a component with temperature in the range 
1.0-2.0 keV. This is consistent with a density gradient 77 
10 — 12. A lower temperature gives too strong Mgxi, and 
Si XIII lines, while a higher temperature gives too strong 
Mgxii and Sixiv lines. However, temperatures above 1.0 
keV strongly overestimate the emission at the Fe L peak at 
0.9 keV, while not having sufficient emission below 0.8 keV. 
This may indicate that there is one component responsible 
for the line emission above 1.0 keV, and another, cooler 
component, responsible for the emission at lower energies. 



For a solar composition the O viii-emission at 0.6- 
0.7 keV is prominent. As this is absent in the observed spec- 
trum, a strong absorption is needed to remove this emission. 
This, however, also removes the N vii-emission at 0.5 keV, 
which is evident in the spectrum. 

We may improve this by adding a second component, 
which gives a best fit for shocks of kT^cv = 1.0 and 
/cTrov = 2.0 keV, shown in Fig [31 We note that while 
the XMM-observations are reasonably well reproduced in 
the range 0.8-1.0 keV, the Chandra observations show a 
dip at 0.9 keV, which means that this model overestimates 
the emission there. Both spectra are poorly reproduced be- 
low 0.7 keV, most likely because of the strong absorption 
needed to remove the O viii-emission. 

We therefore conclude that a composition with less O 
and more N than the solar composition is needed. This 
is in agreement with the results from the optical emis- 
sion, wiiich_2ndicatethe^ of CNO-processed mate- 
rial (jFransson et al.ll2005( ) . In the following we therefore in- 
vestigate the spectra produced by compositions taken from 
the explosion models described above. 

9.2. A CNO enriched hydrogen envelope 

The composition from the outer envelope of the 4H47 
model, where the original envelope composition has been 
modified by GNO-processing in the progenitor (Table [1]), 
gives weaker 0-lines than the solar model, while the N-lines 
are stronger. However, also this composition gives too weak 
Nvii emission at ^ 0.5 keV. These models also show a defi- 
ciency at 0.6-0.7 keV compared to the observed spectrum, 
while the Fe L/Nex peak at 0.9-1.0 keV is too strong. This 
indicates that this composition has too much Ne and too 
little O. At energies above '^1.0 keV the continuum is too 
weak, and there is hardly any line emission. 

The outer envelope of the sl5s7b2f model has more O 
and less Ne than the corresponding zone in the 4H47 model, 
giving stronger O vii-viii lines and a less pronounced Ne x 
peak. Also, the Fexvii lines are stronger, as are the Mgxi- 
XII and Sixiii-xiv lines at 1.2 keV and 1.8 keV. However, 
although this model does a better job at reproducing the 
Fe L emission and the line emission above 1.2 keV, also 
in this model the 0-emission at 0.7-0.8 keV is underes- 
timated, and the Ne-peak is too strong. In this case two 
shocks improve the fit, but cannot reproduce the missing 
0-emission. Fig. [H shows a two-shock model based on the 
model sl517b2f. In this case the temperatures of the two 
shocks are 0.7 and 1.4 keV. 

As a realistic CNO-enriched composition we have also 
tested the SN 1987A-like composition discussed in Sect. [SI 
The fit obtained with this composition is better than with 
the other two CNO-enriched compositions, although the 
Nvii peak at 0.5 keV also in this case is underproduced, 
while the Fexvii emission at 0.8 keV is too strong. 

Adding a second shock improves the fit further by al- 
lowing strong line emission at high energies without over- 
estimating the Fe emission at 0.9 keV. An example, with 
two radiative shocks with Ti = 1.4 keV and T2 — 0.7 keV, 
is shown in Fig. [5l 

In summary, this composition gives an acceptable fit to 
the observations. 
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Fig. 3. Fit to the XMM spectrum of SN 1993J for an adiabatic shock with kT^cv = 1-3 keV (dotted curve) and two 
radiative shock with hTj-^-v = 0.7 keV (upper dashed curve) and fcTj-cv = 0.1 keV (lower dashed curve) for a solar 
composition. The solid curve is the sum of the three components. Note the poor fit at the Nvii peak at ~ 0.5 keV and 
the O VII- VIII emission at ~ 0.7 keV. 



NVII O VII Ne IX-X Mg XI Si XIII S XV ArXVII 
O VIII Fe XIV-XXVI Mg XII Si XIV 




0.5 1 2 5 

channel energy (keV) 



Fig. 4. Fit to the spectrum of SN 1993 J with an adiabatic shock at fcTrov = 1.4 keV (long dashed) and a radiative shock 
at /cTfcv = 0.71 keV (short dashed) for a composition based on the outer, hydrogen-rich envelope of the sl5s7b2f model. 
Also here the fit to the Nvii and Ovii-viii emission at 0.5-0.7 keV is poor. 



9.3. The Helium/Nitrogen zone 



This zone, dominated by CNO burning products, with He 
dominant and N and Ne the most abundant metalsrepro- 
duces the emission around 0.5 keV well. The Fe L and Ne x 
peak at 0.9 keV is, on the other hand, overestimated, and 
there is an obvious lack of emission at 0.6-0.7 keV, due to 
a too low abundance of O. The sl5s7b2f model produces 
an Fe L/Nex peak at slightly lower energy than the 4H47 
model, due to the higher Nex-emission in the 4H47 model. 
Several emission lines are evident above 1 keV, stronger in 
the sl5s7b2f model. 



As with the previous models, adding a second shock 
improves the fit, and the best fit is found for an adiabatic 
shock at fcTj-cv = 2.1 keV and a radiative component at 
fcTrov = 1-0 keV (Fig. O. This gives the best fit to the 
spectrum of all the compositions. 



The highest temperature component is consistent with 
an ejecta density gradient of 77 12 if A* = 4 and V4 = 1 
(Tabled!). 
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Fig. 5. Fit to the XMM spectrum of SN 1993J with a composition resembhng that of the circumstellar ring of SN 1987A, 
but with enhanced metal abundances. The fit consists of an adiabatic shock at kT^-^v — lA keV (dotted) and a radiative 
shock at A;Ti.cv=0.7 keV (dashed). The fit to the Nvii peak at ^ 0.5 keV has improved somewhat compared to Figs[3]|4l 
but the O VII- VIII emission at 0.7 keV is still poorly fit. 



NVII O VII Ne IX-X Mg XI Si XIII S XV ArXVII 
O VIII Fe XIV-XXVI Mg XII Si XIV 
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Fig. 6. Fit to the XMM spectrum of SN 1993 J with an adiabatic shock at fcTicv = 2.1 keV (long dashed) and a radiative 
shock at /cTrov = 1-0 keV (short dashed), for a He/N-dominated composition from the sl5s7b2f model. Note the strong 
improvement in the fit to the Nvii peak at 0.5 keV compared to Figs.[3]-[5l 



9.4. The Carbon/Oxygen zone 



10. Discussion 



A C/O-dominated composition gives spectra that do not in 
any way resemble the observed spectrum. For illustration 
we show one in Fig. [71 Both the N peak at 0.5 keV Fe peak 
at 0.9 keV are completely missing, while the oxygen emis- 
sion at 0.6-0.7 keV and the Mg-S line emission above 1 keV 
are very strong and dominate the spectrum. Changing the 
temperature, or adding more shocks, modifies the strength 
of the peaks, but does not result in an overall improvement. 



In this paper we have investigated the X-ray spectra of 
SN 1993J. By applying a hydrodynamical model for a cool- 
ing shock to the X-ray observations we have discussed the 
fits obtained by various compositions taken from realistic 
ejecta models. 

We find that in the spectra of SN 1993J roughly half the 
contribution comes from an adiabatic shock propagating 
into a CNO-enriched ejecta, with the rest of the contribu- 
tions from slower, radiative shocks. These could be caused 
by instabilities at the reverse shock front or by a clumpy 
CSM modifying the structure of the interaction region. 
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Fig. 7. Fit to the XMM spectrum of SN 1993J with two radiative shocks with fcTrcv = 4.5 kcV (lower dashed curve) and 
fc^rcv — 9.8 (upper dashed curve), for a C/O-dominated composition from the 4H47 model. Note the poor fit at 0.7-1.2 
keV, ruling out this model. 



Comparing ejecta mo dels, we find that the model sl5s7b2f 
(|Wooslev et all 119941) gives a somewhat better fit to the 
line features th an the model 4H47 (jNomoto fc Hashimotol 
I198& .Shigevama et al.lll994D . but the difference is hardly 
significant. The best fits are obtained by a composition re- 
sembling that of the circumstellar ring of SN 1987A. There 
is no need for a contribution from the circumstellar shock, 
as was also expected from the low luminosity estimated in 
Sect, m 

The observations discussed in the present work are close 
to the limit of what can be used for spectroscopic studies. 
If the capabilities of Chandra and XMM had been present 
at the early phases of this supernova, a much more detailed 
analysis could have been performed of the X-ray spectra. 
The next X-ray bright supernova within ~ 10 Mpc will al- 
low a much more detailed spectroscopic studies of the early 
spectra. With future missions, like ESA's, NASA's and 
JAXA's proposed International X-ray Observatory (IXO), 
these capabilities will be vastly improved. This will allow 
high S/N spectra to be obtained of SN at much longer dis- 
tances, and will also allow detailed spectroscopic studies 
to be performed also at late times. This will in turn lead 
to a much more reliable determination of the abundances 
and temperature in the interaction region. We hope that 
the analysis in this paper will give some guidance to future 
X-ray observations of SNe. 
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Table 1. Fractional abundances by number and molecular weights for the compositions discussed here. 
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3 










5.15 X 10- 


-4 





6.29 


X 


10" 


4 


3.38 


X 


10" 


5 


5.08 


X 


10" 


-5 


4.80 


X 


10" 


-1 


6.99 


X 


10" 


-4 


6.49 


X 


10" 


5 


7.43 


X 


10" 


-1 


4.69 X 10" 


-4 


Ne 


8.89 


X 


10" 


5 


7.86 


X 


10" 


4 


1.77 


X 


10" 


-3 


1.23 


X 


10" 


-2 


1.19 


X 


10" 


-4 


3.56 


X 


10" 


4 


1.93 


X 


10" 


-2 


8.89 X 10" 


-5 


Mg 


3.62 


X 


10" 


-5 


3.57 


X 


10" 


-5 


8.03 


X 


10" 


-5 


3.21 


X 


10" 


-3 


3.92 


X 


10" 


-5 


1.17 


X 


10" 


4 


7.69 


X 


10" 


-4 


3.63 X 10" 


-5 


Si 


3.46 


X 


10" 


-5 


4.42 


X 


10" 


-5 


9.94 


X 


10" 


-5 


3.25 


X 


10" 


-4 


3.73 


X 


10" 


-5 


1.11 


X 


10" 


4 


4.16 


X 


10" 


-4 


3.45 X 10" 


-5 


S 


1.69 


X 


10" 


-5 


2.13 


X 


10" 


-5 


4.78 


X 


10" 


-5 


1.56 


X 


10" 


-4 


1.81 


X 


10" 


-5 


5.40 


X 


10" 


6 


2.02 


X 


10" 


-4 


1.69 X 10" 


-5 


Ar 


3.62 


X 


10" 


6 


9.30 


X 


10" 


6 


2.08 


X 


10" 


-5 


6.13 


X 


10" 


-5 


3.63 


X 


10" 


-6 


1.08 


X 


10" 


5 


4.05 


X 


10" 


-5 


3.62 X 10" 


-6 


Ca 


2.13 


X 


10" 


6 


2.63 


X 


10" 


6 


5.93 


X 


10" 


-6 


1.94 


X 


10" 


-5 


2.23 


X 


10" 


-6 


6.66 


X 


10" 


6 


2.49 


X 


10" 


-5 


2.13 X 10" 


-6 


Fe 


3.08 


X 


10" 


5 


5.30 


X 


10" 




1.19 


X 


10" 


-4 


3.90 


X 


10" 


-4 


3.28 


X 


10" 


-5 


9.79 


X 


10" 


5 


3.66 


X 


10" 


-4 


3.08 X 10" 


-5 


Ni 


1.62 


X 


10" 


6 


2.78 


X 


10" 


6 


2.22 


X 


10" 


-6 


2.05 


X 


10" 


-5 


1.72 


X 


10" 


-6 


5.14 


X 


10" 


6 


1.92 


X 


10" 


-5 


1.62 X 10" 


-6 


HA 


1.29 








1.81 








4.07 








11.37 








1.36 








4.06 








15.20 








1.63 






0.61 








0.80 








1.34 








1.70 








0.64 








1.34 








1.77 








0.74 





Table 2. Reverse shock temperature, fcTrev and limiting temperature, fcTrad (both in keV) for a radiative shock, for 
SN 1993J at 2600 days for M = 4.0 x lO-^Mg yr"!. 



n 


7 


12 


20 


Composition 


kTrad 




kTrad 








Solar 


0.51 


4.77 


0.63 


1.19 


0.68 


0.37 


H/Hc (4H47) 


0.69 


6.30 


0.82 


1.57 


0.88 


0.48 


Hc/N (4H47) 


0.93 


10.5 


1.04 


2.60 


1.09 


0.81 


C/O (4H47) 


2.02 


13.5 


2.43 


3.39 


2.61 


1.04 


H/He (sl5s7b2f) 


0.53 


5.01 


0.66 


1.25 


0.71 


0.39 


He/N (sl5s7b2f) 


0.86 


10.4 


0.99 


2.60 


1.04 


0.80 


C/O (sl5s7b2f) 


2.32 


13.9 


2.80 


3.46 


3.00 


1.07 


SN 1987A-like 


0.73 


5.70 


0.87 


1.40 


0.92 


0.44 



Table 3. Fitting parameters for different components for the SN 1993 J XMM-Newton data. 



Parameters 


Models 


Solar 


HHe-Zonc 

(sl5s7) 


HcN-Zonc 

(sl5s7) 


SN87A-Uke 


x~/<H' 


2.i:i 


i.98 


i.2() 


i.()r 


Rad. Comp ^cmp (kcV) 

Norm (10^^) 


0.68 ±0.08 
5.79 ±0.11 


0.71 ± 0.08 
6.53 ± 0.83 


1.02 ±0.41 
6.20 ±0.23 


0.75 ±0.11 
5.10 ±0.73 


Ad. Comp T^^^P ('^'^y^ 

Norm (lO'^) 


1.33 ±0.24 
2.08 ± 0.04 


1.37 ±0.03 
1.81 ±0.05 


2.07 ±0.16 
0.28 ± 0.02 


1.35 ± 0.03 
1.58 ± 0.03 


Unabsorbed flux * Rad. Comp 
10"^^ erg s"^ Ad. Comp 

Total flux 


1.12 ±0.02 
1.72 ±0.03 
2.84 ± 0.05 


1.24 ±0.15 
1.59 ±0.04 
2.84 ±0.19 


1.64 ±0.6 
1.45 ±0.10 
3.09 ±0.16 


1.14 ±0.16 
1.73 ± 0.03 
2.87 ±0.19 


nH(cool shell) in 10''" cm"'' 


13.2 ± 1.40 


11.8 ±0.96 


0.07 ± 0.07 


6.5 ±1.1 


Flux (Reverse Shock) ** Rad. Comp 
10"^^ erg s"^ Ad. Comp 

Total flux 


2.49 ±0.05 
2.60 ± 0.05 
5.09 ±0.10 


2.46 ± 0.31 
2.92 ±0.08 
5.38 ± 0.39 


2.73 ±0.10 
1.89 ±0.14 
4.62 ±0.24 


2.18 ±0.31 
2.46 ± 0.05 
4.66 ± 0.36 



* Corrected for the Galactic extinction column density, nn~ 4.17 x 10^" cm 
** The total flux from the reverse shock before absorption by the cool shell 



